Journal of Modern Mathematics Frontier Vol. 1 Iss. 4, December 2012 



A Theoretical Model for an Immobilized 
Oxidase Enzyme Electrode in the Presence of 
Two Oxidants 

Anitha Shanmugarajan, Rajendran Lakshmanan 

Department of Mathematics, The Madura College, Madurai-625011, Tamilnadu, India 
raj_sms@rediffmail.com 



Abstract 

A theoretical model of an immobilized oxidase enzyme 
electrode is discussed. The model is based on three 
reaction/diffusion equations containing a non-linear reaction 
term under the steady state conditions. An analytical 
expression pertaining to concentrations of the 
immobilization of three enzyme substrates are obtained. 
Furthermore, we employ Homotopy perturbation method 
(HPM) to solve the system of non-linear reaction/dif fusion 
equations. Simple and an approximate polynomial 
expression of concentration of substrate, oxygen and 
oxidized mediator and current were obtained in terms of the 
Thiele moduli and the small values of parameters B s 
(normalized surface concentration of substrate), B a 
(normalized surface concentration of oxygen) and B m 
(normalized surface concentration of oxidized mediator). 
The minimum value of the dimensionless concentration of 
F s , F and F m in terms of parameter is also reported. 

Furthermore, in this work the numerical simulation of the 
problem is also reported using Matlab program. An 
agreement between analytical expressions and numerical 
results is observed. 
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Introduction 

The concept of enzyme electrode was first introduced 
by Clark and Lyons [1] over fifty years ago. Various 
kinds of biosensors based on electrochemical enzyme 
electrodes have been studied and a majority of these 
devices operates in an amperometric mode [2]. 
Schulmeister [3, 4] described a model for a modified 
amperometric enzyme electrode with a multilayer 
enzyme membrane. This model is employed to 
describe the behavior of the electrode in enzyme 
kinetics. The kinetics is described by a parabolic 
differential equation with linear inhomogeneities. 



Leypoldt and Gough [5] have proposed a 
mathematical model for glucose oxidase immobilized 
enzyme electrode for various diffusion coefficients, 
and kinetic parameters corresponding to an enzymatic 
reaction. Bergel and Comtat [6] studied the transient 
response of a mediated amperometric enzyme 
electrode using an implicit finite difference method. A 
model for two-substrate enzyme electrode was 
described by Gooding and Hall [7]. 

The importance of enzyme electrode and the various 
approaches to enzyme electrochemistry has been 
discussed. The current response has been also 
observed experimentally for mediated enzyme 
electrodes employing glucose oxidase and other 
oxidases [8-10]. In recent years homogeneous 
membranes are widely used as carriers for 
immobilization of enzymes [11]. They have been 
utilized in biomaterials, bio-separators and biosensors 
[12]. The membranes provide an ideal support for the 
immobilization of the biocatalyst. Substrate partition 
at the membrane/fluid inter-phase can be used to 
improve the selectivity of the catalytic reaction 
towards the desired products [13]. 

Recently, Martens and Hall [14] have described a 
model for three-substrate homogeneous enzyme 
electrode and employ the model to investigate the 
influence of oxygen on the current response of a 
mediated oxidase electrode. To our knowledge, no 
general analytical expressions that describe the 
concentration of substrate, oxygen and oxidized 
mediator for various values of the Thiele moduli 
S 2 , 0l and 2 and the normalized parameters 
B s , B Q and B m have been reported. However, in 
general, analytical solutions of non-linear differential 
equations are more important and very useful than 
purely numerical solutions, as they are amenable to 
various kinds of data manipulation and analysis. For 
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this reason, we have derived that analytical 
expressions corresponding to the concentrations of 
substrate, oxygen and oxidized mediator in an oxidase 
enzyme electrode using Homotopy perturbation 
method. This method is reliable and highly accurate in 
handling non-linear problems. 

Mathematical Formulation of the Boundary 
Value Problem 

The details of the model adopted have been fully 
described in Martens and Hall [14] and so we only 
present a brief summary here. Fig. 1 represents the 
general kinetic reaction scheme of enzyme-membrane 
electrode geometry [15]. 

Far wall Membrane Electrode Solution 




- 



i d 

FIG. 1 SCHEMATIC MODEL OF AN ENZYME-MEMBRANE 
ELECTRODES [12] 

Eqs. 1-3 represent the general reaction scheme that 
occurs often at an electrode immobilized with an 
enzyme oxidase in the presence of two oxidants: 



E ox +S^-+ES^-+E red +P 



E red +0, 



E red +Med ox 



->E„ +H,0 



2 W 2 



^E ox +Med red 



(1) 

(2) 
(3) 



We assume that the concentrations of all reactants and 
enzyme intermediates remain constant for all time. 
Also the concentration of total active enzyme [E t ] and 
the reactants in the bulk electrode remain constant. We 
can consider that the diffusion of the reactants can be 
described by Fick's second law and the enzymes are 
assumed to be uniformly dispersed throughout the 
matrix. The enzyme activity is not a function of 
position. The reaction/dif fusion equations 
corresponding to the concentrations of substrate, 
oxidized mediator and oxygen within a matrix can be 
expressed as [14] 
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where [S], [Med | and [O2] are the concentrations 

of substrate, oxidized mediator and oxygen and y is 
distance from the electrode. D s , D m and D are the 
diffusion coefficients of substrate, oxidized mediator 
and oxygen. [E t ] Denotes the concentration of the 

total active enzyme in the matrix and /? s = (k 2 + k i )/k l , 
p = k 3 /k 4 & P m = k 3 /k 5 represent reaction constants 
of substrate, oxygen and oxidized mediator 
respectively. The boundary conditions are given by 

When y = d , 



[0 2 ] = [0 2 ] b =K [0 2 ] a , [S] = [s] b =tf s [s] b , 



[Med ox ] = [Med ox ] b =^ m [Med ox ] a0 



(5a) 



when v = , 

[Med ox ] = [Med ox ] b , d[o 2 ]/dy =d[s]/dy =0 (5b) 

Here [S]„, [OJ b , [MedJ b & [S]., [0 2 ]., [Med ox ] 00 

are the bulk concentrations of substrate, oxygen and 
oxidized mediator species. K a , K s and K m are the 
partition coefficients for oxygen, substrate and 
oxidized mediator respectively. Eq. 4 can be 
represented in the normalized form using the 
following normalized parameters. 
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[Med ox ] b 



fin 



(6) 
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' ' ' y o ' 



D ot0 2 ] b 



d k 3 [E t ] 
fl m [ Med ox] b 



where F s , F , and F m represent the normalized 

concentrations of substrate, oxygen and oxidized 
mediator and Z? s , B a , and B m are the corresponding 

normalized surface concentrations. The surface 
concentration is the ratio of the bulk concentration and 
the reaction constants. S ,0 O , and @ m denote the 
Thiele moduli of substrate, oxygen and oxidized 
mediator, respectively. Thiele modulus @ 2 represents 
the ratio of the characteristic time of the enzymatic 
reaction to that of substrate diffusion, d is the thickness 
of the enzyme layer. The coupled three non-linear 
reaction/diffusion equations in normalized form are 



is a reliable and efficient method to solve the non- 
linear problems [18-24]. The basic principle of this 
method is described in [21]. Detailed derivations of 
dimensionless concentration profiles are described in 
Appendix A. As a result, we have obtained the 
analytical expression corresponding to the 
concentrations of substrate, oxygen and oxidized 
mediator as follows: 



F s (x)=<Z> 2 (l-x 2 ) 



2 + - 



B n + B„ 



B, 



(x)=1 _ ^o(l-x 2 ) 
M) (B 0+ B m ) 

(x-x 2 ) 



2 B 



2 + 



B +B„ 
2 



B t 



B n + B n 



B. 



(10) 



(11) 



(12) 



From Eq. 9, we can obtain the current as follows: 



<t> 2 B 



(b +bJ 



2 + - 



B n +B„ 



B. 



(13) 
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F B + F m B m F S B S j 



Eqs. 10-12 represent the new analytical expression of 

(7a) 

concentrations of substrate, oxygen and oxidized 
mediator. 
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1 d rF m F B +F m B m 



<t>l dx 2 



F B +F m B m F S B S j 



(7c) 



The boundary conditions become 
F . = F = F m = 1 at X = 1 



dF. 6F. „ 
: 1, —5- = — ^ = 0atX = 
dx dx 



The normalized current J ox is given by, 



(8a) 
(8b) 



dx 



(9) 



J x=0 



Analytical Expressions of Concentrations of 
Substrate, Oxygen and Oxidized Mediator 
under Steady-state Condition 

Recently, Ji-Huan He solved many non-linear 
diffusion equations using Homotopy perturbation 
method [16, 17]. The Homotopy perturbation method 



Numerical Simulation 

The three reaction/diffusion equations (Eqs. 7a-7c) 
corresponding to the boundary conditions (Eqs. 8a and 
8b) were solved by numerical methods. We have used 
pdex4 to solve these equations (Pdex4 in MATLAB is a 
function to solve the initial-boundary value problems 
of differential equations). The numerical solution is 
compared with our analytical results and is shown in 
Fig. 5. A satisfactory agreement is noticed for various 
values of the Thiele moduli 2 , 2 , and 0^ and 
possible small values of parameters B s , B and B m . 
Tables 1-3 illustrate the comparison of analytical result 
obtained in this work with the numerical result. The 
relative difference between numerical and analytical 
results is within 3% for all possible values of 
parameter considered in this work. 

Discussion 

The kinetic response of a biosensor depends on the 
concentrations of substrate, oxygen, and mediator. The 
concentrations of substrate, oxygen and mediator 
depends on the following four factors 
2 , fi s , B a and B m . The variation in the Thiele 
modulus can be achieved by varying either the 
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thickness of the enzyme layer or the amount of 
enzyme immobilized in the matrix. The Thiele 
modulus is indicative of the competition between the 
diffusion and reaction in the enzyme layer. When <t> 2 
is small, the kinetics dominate and the uptake of 
oxygen, mediator and substrate in the enzyme matrix 
is kinetically controlled. Under these conditions, the 
substrate concentration profile across the membrane is 
essentially uniform. The overall kinetics are governed 
by the total amount of active enzyme [E t \. Diffusion 
limitations are the principal determining factor when 
Thiele modulus is large. 

A series of normalized concentration profiles for 
substrate F s , oxygen F a , and mediator F m for 

various values of diffusion coefficient of substrate D s , 
oxygen D D , and mediator D m respectively are 
presented in Fig. 2. From Figs, 2a and 2b, it is 
concluded that the concentration profiles of substrate 
and oxygen increases and reaches the maximum value 
when the diffusion coefficients of substrate and 
oxygen increases. The concentration of oxidized 
mediator starts from the value of unity, goes through a 
minimum and then increases back to unity as the 
diffusion coefficient of mediator is increased (Fig. 2c). 






FIG. 2 NORMALIZED CONCENTRATION PROFILES FOR (a) THE 
SUBSTRATE F s , (b) THE OXYGEN F Q , AND (c) THE OXIDIZED 
MEDIATOR F m . THE CURVES ARE PLOTTED USING EQS. 10 - 
12 FOR VARIOUS VALUES OF THE DIFFUSION COEFFICIENTS 
OF SUBSTRATE, OXYGEN, AND OXIDIZED MEDIATOR AND 

SOME FIXED VALUES OF PARAMETERS B s = 0.5 , B Q = 0.01 , 
B m =0.1 , d = 0.01 CM, k 3 = 0.01 LMOL-'S" 1 , [£,] = 0.005 



MMOL/L, AND [S], 



[O, 



[M ] =0.5 MMOL/L. SOLID 
red b 



J b 2 

LINES REPRESENT THE ANALYTICAL SOLUTION PRESENTED 
IN THIS WORK AND DOTTED LINES THE NUMERICAL 
SIMULATION 



The normalized concentration profiles of substrate F s , 
oxygen F D , and mediator F m for various values of 
the normalized reaction constants 

A = (*2 + k 3 ' A> = k 3 l k A & Pm = k 3 l k S shown 

in Fig. 3. From Fig. 3, all curves increase as the reaction 
constants increase while keeping the other parameter 
fixed. From this figure, it is evident that the 
concentration profiles attains the steady state when the 
value of reaction constants > 100 . Fig. 4 represent the 
normalized concentration profiles for various values 
of the Thiele moduli c5 s 2 , cP 2 , and @^ and some fixed 
values of the normalized parameters fi s , B and B m . 
From Fig. 4a, it is evident that the concentration of 
substrate decreases when Thiele modulus increases. 




FIG. 3 A PLOT OF NORMALIZED CONCENTRATION PROFILES 
FOR (a) THE SUBSTRATE F s , (b) THE OXYGEN F , AND (c) THE 

OXIDIZED MEDIATOR F m . THE CURVES ARE PLOTTED USING 
EQS. 10 - 12 FOR VARIOUS VALUES OF THE REACTION 
CONSTANTS AND SOME FIXED VALUES OF 

PARAMETERS &} = <Z>q = <Z>m = 5 , [S] b = 0.01 MMOL/L, 

[0 2 ], =0.001 MMOL/L, AND [M rgd ], = 0.05 MMOL/L. SOLID 

LINES REPRESENT THE ANALYTICAL SOLUTION PRESENTED 
IN THIS WORK AND DOTTED LINES THE NUMERICAL 
SIMULATION 

It is obvious that the substrate concentration 
F s reaches the maximum value 1, at membrane 
/solution interface (x = l). From Fig. Ab, it is evident 
that the value of oxygen concentration F « 1 for all 

small values of cP 2 and large values of k 2 and k 3 . 
Furthermore, there is a simultaneous decrease in the 
values of concentration of oxygen when the Thiele 
modulus increases. The concentrations of F s and F a 

have minimum values 1 - <t>l [2 + 2/(B + B m ) + 2/B s ] _1 
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FIG. 4 NORMALIZED CONCENTRATION PROFILES FOR (a) THE 
SUBSTRATE F s , (6) THE OXYGEN F , AND (c) THE OXIDIZED 

MEDIATOR F m IN THE POLYMER MATRIX FOR VARIOUS 
VALUES OF THE THIELE MODULI <Z>;, <Z> 2 and <Z>* WHEN 
B s = 0.5, S D = 0.05 AND B m = 0.1 . SOLID LINES REPRESENT 
THE ANALYTICAL SOLUTION PRESENTED IN THIS WORK 
AND DOTTED LINES THE NUMERICAL SIMULATION 

and l-^Wk +B m) +2 / B sV '( B o +«m) at 
x = 0. Upon careful evaluation of Fig. 4c, it is evident 
that there is a simultaneous decrease in the values of 
concentration of oxidized mediator when the Thiele 
modulus <X>m or thickness of the enzyme layer 
increases. The concentration F m decreases initially 
and attains its minimum value 
l-0,;B m [2 + 2/(B o +B m )+2/Bj~ l /{B +B m ) at X = x h 
(middle of the membrane) and then increases. Upon 
comparison of analytical with the numerical, it is 
evident that both results are identical. Figure 5 is the 




20 40 eo eo 100 



FIG. 5 THE NORMALIZED CURRENT J m VERSUS THE THIELE 
MODULUS 0l . THE CURRENT WERE COMPUTED USING EQ. 
13 FOR VARIOUS VALUES OF THE NORMALIZED 
PARAMETERS (a) B s =1, B„ = 0.01, B m =10 , (b) 

B s = 10, B a = 1, B m = 0.5 , (c) B s = 0.1, B = 1, B m = 10 , AND (d) 
B s =l, B o =10, B m =0.1. 



normalized current response / for various values of 

the normalized parameters using Eq. 13. From Fig. 5, it 
is evident that the value of the current increases when 
0^ increases. Fig. 6 indicates the profiles of 




FIG. 6 PROFILES OF \ VERSES SURFACE CONCENTRATION 

<Z>; 

OF (a) SUBSTRATE fi s , (b) OXIDIZED MEDIATOR B m , AND (c) 
OXYGEN B IN THE POLYMER MATRIX. THE CURVES ARE 
COMPUTED USING EQ. 13 FOR VARIOUS VALUES OF THE 
NORMALIZED PARAMETERS B m , B and B s 

J ox / <Z>^ for various values of the normalized 
parameters B s , B D and B m using Eq. 13. From Fig. 6a, 
it is clear that the value of the current increases when 
/? m (= k 2 /k 5 ) decreases. In addition, we noticed that 

the normalized current reaches the steady-state value 
when B m <0.1 and Z? s >4. From Fig. 6b, it is obvious 

that the value of the current decreases when 
fi {=k 3 /k 4 ) increases. Furthermore, the value of the 
current reaches the steady-state value when 
B D < 0.01 and B m > 4 . From Fig. 6c, it is observed that 
the value of the current increases when 
Ps ( = (^2 + ^3 )Ai ) decreases and reaches the steady- 
state value when Z? s < 0. 1 and B > 4 . 

Tables 1-3 represent the comparison of analytical 
result with the numerical result of normalized 
concentration profiles for various values of the Thiele 
modulus. From the table 1, it is inferred that the 
concentration F s decreases when c£> s 2 increases or the 
thickness of the electrode increases. From the table 2, it 
is evident that the normalized substrate concentration 
F w 1 when the value of < 1 . The comparison 
reveals that the relative difference between the 
analytical results and the numerical solution does not 
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exceed 3% for all values of parameters considered in 
the simulation. 



Here p e [0, 1] is an embedded parameter. Applying Eq. 
A2 for the Eqs. 7a-7c, we get 



Conclusions 

We have analyzed the theoretical model for an 
immobilized oxidase enzyme electrode, previously 
described in Martens and Hall [14]. The system of non- 
linear steady-state reaction/diffusion equations of the 
model has been solved analytically. Approximate 
analytical expressions pertaining to the concentrations 
of species and current in three-substrate enzyme 
electrodes for all values of the Thiele moduli are 
obtained using the Homotopy perturbation method. 
Our results agree well with the numerical solution. 
This analytical result is useful for improve the sensors 
design. The extension of the procedure to the model 
with reduced mediator in the bulk solution seems 
possible. The theoretical results obtained can be used 
for the optimization of the thickness of the enzyme 
layer (or Thiele modulus) which produce significant 
change in both the magnitude of the current and the 
general behavior of the system. 
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Appendix A 
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where linear operator L(F S ) = d 2 F s / dx 2 , 
L(F ) = d 2 F /dx 2 , and L(F m ) = d 2 F m I dx 2 and 

non-linear operator 



N(F S ) = 



1 + 1 + - 1 



F B o +F m B m f Ay 



Analytical solution of the system of non-linear Eqs. 
7a-7c using Homotopy perturbation method 

In this appendix, we derive the solution of non-linear 
reaction/diffusion equations Eqs. 7a-7c using He's 
Homotopy perturbation method. To illustrate the 
basicconcepts of this method (HPM), we consider the 
following non-linear differential equation 



N(F ) = 



and 



N(F m )= 



1 + + 

V F B o +F m B m 



Fo B o 



F o B o +F m B m 
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F m B m 



F o B o +F m B m 



L(u) + N(u)-f(r) = , 



n 



(Al) 



where L is a linear operator, N is a nonlinear operator, 
u is an unknown function, and fir) is a given 
continuous function. We construct a homotopy which 
satisfies 



(1 - p)[L(u)]+ p[L(u) + N(u) -f(r)] = 



(A2) 



Supposing the approximate solutions of Eqs. A3 - A5 
have the form 



F s = F sfi + P F s,l + P F s,2 + ■■ 
F o = F o,Q + P F o,l + P 2 F o,2 +• 



= F m,0 + P F m,l + P F m,2+- 



(A6) 
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Substituting Eq. A6 into Eqs. A3 - A5 (respectively), 
and equate the terms with the identical powers of p, 
we obtain 



p ■ 



dx 1 



(A7) 
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P ■ t~ n \ = (A8) 

dx 2 1 + hVs 
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We can obtain the coefficient of p for F as in the 
Eq. A9. Now substituting the third equation of A6 in 
A5, we can obtain 



d 2 F, 



m,l 



(A12) 



1 



dx 2 l + l/{F ofl B + F mfl B s )+F sfl B s 



F m,0 B m 



(A13) 



Similarly we can find the coefficient of p for F m as 
in the Eq. A9. The initial approximations (boundary 
conditions Eqs. 8a and 8b) are as follows: 

F s0 (x = l) = l; F o() (x = l) = l; F m0 (x = l) = l (A14) 
F m Q (x = 0) = 1; dF Q fi (x = Q)jdx = dF Q (x = 0)/dx = 

(A15) 

and 

F si (x = l) = 0;F oi (x = l) = 0; F mi (x = l)=0; 

i =1,2,3,-. •• (A16) 



F m . (x = 0) = ; dF 0ji (x = 0)/dx =dF Sji (x = 0)/dx = ; 

f=l,2,3,... (A17) 

Solving the Eqs. A7, A10, and A12 and using the 
boundary conditions Eqs. A14 and A15, we get 



0,0 



and 



(A18) 
(A19) 

(A20) 



Substituting the above values of F s0 , F o0 and F mfi in 

the Eq. A8, Eq. All and Eq. A13 and solving these 
equations, using the boundary conditions Eqs. A16 
and A 17, we can obtain the following results. 

.2 



F s,l = 
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04 [2 + 2/(B 0+J B m )+2/ J B s ] 
and 



B„+B„ 



0- 



m,l 



[2 + 2/(B 0+ B m )+2/B s ] 



B +B n 



(x 2 -l) (A22) 



(x 2 -x) (A23) 



Adding Eq. A18 and A21, we get Eq. 10 in the text. 
Similarly we can get Eqs. 11 and 12 in the text. In Eq. 
A6, we can consider only two terms (two iterative 
steps). To find the next iterative, we can solve the Eq. 
A9, after substituting the Eqs. A21-A23. Accuracy of 
the solutions can be improved by considering next 
terms ( ie by solving the Eq. A9). 

Appendix B 

Nomenclature 
Symbols 

[E t ] -Total active enzyme concentration in the 

matrix (mmol/L) 
[E ox ] -Enzyme concentration of the oxidized 

mediator (mmol/L) 
[ES] -Enzyme concentration of the substrate (mmol/L) 
[E red ] -Reduced enzyme concentration (mmol/L) 
[0 2 ] -Concentration of oxygen at any position 
in the enzyme layer (mmol/L) 
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[0 2 ] b -Oxygen concentration in the bulk electrolyte 
(mmol/L) 

[0 2 ] K -Oxygen concentration in the bulk solution 
(mmol/L) 

[S] -Concentration of substrate at any position in the 
enzyme layer (mmol/L) 

[S] b -Substrate concentration in the bulk electrolyte 
(mmol/L) 

[S]^ -Substrate concentration in the bulk solution 
(mmol/L) 

[Med ox ] -Concentration of oxidised mediator at any 

position in the enzyme layer (mmol/L) 
[Med ox ] b -Oxidised mediator concentration in the bulk 

electrolyte (mmol/L) 
[Med ox ] x -Oxidised mediator concentration in the 

bulk solution (mmol/L) 

D -Diffusion coefficient of oxygen (cm 2 s- 1 ) 

D s -Diffusion coefficient of substrate (cm 2 s- 1 ) 

D m - Diffusion coefficient of mediator (cm 2 s _1 ) 

d -Enzyme layer thickness (cm) 

y -Distance from the electrode (cm) 

k lr k 4 , k 5 - Rate constants ( L mob 1 s _1 ) 

k 2 , k 3 -Rate constants ( s 4 ) 

K Q -Partitioning coefficient for oxygen (none) 

K s - Partitioning coefficient for substrate (none) 

K m -Partitioning coefficient for mediator (none) 

B Q -Normalized surface concentration of oxygen 
(none) 

B s - Normalized surface concentration of the 
substrate (none) 

B m - Normalized surface concentration of mediator 
(none) 

F - Normalized oxygen concentration (none) 
F s - Normalized substrate concentration (none) 
F m - Normalized mediator concentration (none) 



x - Normalized distance from the electrode (none) 
J m - Normalized current (none) 

Greek symbols 

<j>l - Thiele modulus for the oxygen (Normalized) 

(frl - Thiele modulus for the substrate (Normalized) 

(j) 2 m - Thiele modulus for the mediator (Normalized) 
Subscripts 



o - Oxygen 

s - Substrate 

m - Mediator 

ox - Oxidized species 

red - Reduced species 

t - Total 

oo - Bulk solution 
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TABLE. 1 COMPARISON OF NORMALIZED SUBSTRATE CONCENTRATION F & WITH THE NUMERICAL RESULT FOR VARIOUS 
VALUES OF <Z> 2 AND SOME FIXED VALUES OF NORMALIZED PARAMETERS B s = 0.04, B Q = 2.5 AND B m = 0.2 . 
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TABLE. 2 COMPARISON OF NORMALIZED OXYGEN CONCENTRATION F Q WITH THE NUMERICAL RESULTS FOR VARIOUS 
VALUES OF 1 AND SOME FIXED VALUES OF NORMALIZED PARAMETERS B s = 0.04, B Q = 2.5 AND B m = 0.2 . 
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TABLE. 3 COMPARISON OF NORMALIZED MEDIATOR CONCENTRATION F m WITH THE NUMERICAL RESULT FOR VARIOUS 
VALUES OF <Z> m AND SOME FIXED VALUES OF NORMALIZED PARAMETERS B s = 0.07, B Q = 2.5 AND B m = 0.2 . 
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TABLE. 4 NUMERICAL VALUES OF THE PARAMETERS USED IN THIS WORK. THE FIXED VALUES OF THE PARAMETERS ARE 
[S] = 5mmol/L,[0 2 ] b = 1 .2 mmol/L , [Med ] = 1 mmol/L, AND[£ t ] = 0.26 mmol/L, d = 1 x 10~ 3 cm, k y = 14,000 Lmol" 1 * -1 , 

k 2 = Lmof 1 .?" 1 ,^ = 1,000 Lmol" 1 ^ -1 , = 2 x 10 6 Lmof \v~ 1 , ^ =2x 10 5 Lmof V?" 1 , P s = 0.0714 MOL/L, P Q = 0.0005 

MOL/L, fi m = 0.005 MOL/L. THESE PARAMETERS ARE USED IN MARTENS AND HALL [14], 
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